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ABSTRACT 



Context. The Helios measurements of the angular momentum flux L of the fast solar wind lead to a tendency for the fluxes associated 
with individual ion angular momenta of protons and alpha particles, L,, and L„, to be negative (i.e., in the sense of counter-rotation 
with the Sun). However, the opposite holds for the slow wind, and the overall particle contribution Lp = L,, + La tends to exceed the 
magnetic contribution Lm- These two aspects are at variance with previous models. 

Aims. We examine whether introducing realistic ion temperature anisotropies can resolve this discrepancy. 

Methods. From a general set of multifluid transport equations with gyrotropic species pressure tensors, we derive the equations 
governing both the meridional and azimuthal dynamics of outflows from magnetized, rotating stars. The equations are not restricted 
to radial flows in the equatorial plane but valid for general axisymmetric winds that include two major ion species. The azimuthal 
dynamics are examined in detail, using the empirical meridional flow profiles for the solar wind, constructed mainly according to 
measurements made in situ. 

Results. The angular momentum flux L is determined by the requirement that the solution to the total angular momentum conservation 
law is unique and smooth in the vicinity of the Alfven point, defined as where the combined Alfvenic Mach number Mr = 1. Mr has 
to consider the contributions from both protons and alpha particles. Introducing realistic ion temperature anisotropies may introduce 
a change of up to 10% in L and up to ~ 1.8 km s ' in azimuthal speeds of individual ions between 0.3 and 1 AU, compared with the 
isotropic case. The latter has strong consequences on the relative importance of Lp and Lm in the angular momentum budget. 
Conclusions. However, introducing ion temperature anisotropies cannot resolve the discrepancy between in situ measurements and 
model computations. For the fast-wind solutions, while in extreme cases Lp may become negative, L,, never does. On the other hand, 
for the slow solar wind solutions examined, Lp never exceeds L^, even though may be less than the individual ion contribution, 
since L,, and L„ always have opposite signs for the slow and fast wind alike. 

Key words. Sun: rotation - Sun: magnetic fields - solar wind - Stars: rotation - Stars: winds, outflows 



1. Introduction 



The angular momentum loss of a rotating star due to its 
outflow influences the rotational evolution of the star con- 
siderably, and is therefore of astrop h ysical significance in 
general (see e.g., Weber & Davisl 1967: Belc her & MacGregoJ 
ll976HMestel & Spruillll987tlBouvieret al.iil997l) . However, di- 
rect tests of in si tu measuremen t s aga inst theories such as 
those presented bv I Weber & Davisl d 19671) are only possible for 
the present Sun. A substantial number of studies have been 
conducted and we re compiled in the comprehensive paper by 
iPizzo et all (Il983b . who themselves paid special attention to the 
Helios measurements of specific ang ular momentum fluxes. Th e 
measurements, further analyzed by i Marsch & Richteij ( 11984 . 
are unique in that they allow the individual ion contribution 
from protons Lp and alpha particles L,, to the solar angular- 
momentum loss rate per steradian L to be examined. For in- 
stance, despite the significant scatter, the data exhibit a distinct 
trend for Lp to be positive (negative) for solar winds with pro- 
ton speeds Vp below (above) 400 km s"'. A similar trend for 
La is also found on average. The magnetic contribution Lm, 
on the other hand, is remarkably constant. A mean value of 
Lm = 1.6 X lO^** dyne cm sr ' can be quoted for the so- 
lar winds of all flow speeds and throughout the region from 
0.3 to 1 AU. For comparison, the mean values of angular mo- 



mentum fluxes carried by ion flows in the slow solar wind are 
L„ = 19.6 and L„ = 1.3 xlO^'' dyne cm sr ' (see Table II 
of iPizzo et al.lll983h . The overall particle contribution to L is 



then Lp = Lp + La = 20.9 x lO''' dyne cm sr , which tends 
to be larger than Lm- It is noteworthy that a more recent study 
by IScherer et alj ( 1200 lb showed how examining the long-term 
variation of the non-radial components of the solar wind veloc- 
ity and the corresponding angular momentum fluxes can help us 
understand the heliospheric magnetic field better. 

Alpha particles should be placed on the same footing as pro- 
tons from the perspective of solar wind modeling, given their 
non-negligible abundance and the fact that there tends to exist 
a substantial differential speed v„p = |Vo,/,|sign(|v„| - |Vp|). As 
shown by the Helios measurements, a Vap amounting to up to 
20% - 30% of the local proton speed may oc cur in both the 
fast and slow solar winds (iMarsch et al.lll982allbl) . with the lat- 
ter being exemplified by an event that took place on day 117 
of 1978, when a positive V a„ ~ 100 km s ' was found at 
0.3 AU (IMarsch et al.lll981h . That on the average v„p =s in 
the slow wind simply reflects that the events wit h positive and 
negativ e Vap occur with nearly equal frequency dMarsch et al.l 
Il982al) . As for the alpha abundance relative to protons, a value 
of 4.6% (0. 4%- 10%) is well-esta bHshed for the fast (slow) solar 
wind (e.g.. lMcComas et"aril200Cll) . Therefore alpha particles can 
play an important role as far as the energy and linear momen- 
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turn balance of the solar wind are concerned. When it comes to 
the problem of angular momentum transport, it was shown that 
in interplanetary space not only the angular momentum flux car- 
ried by the alpha particles L„ but also that convected by the pro- 
tons L „ are determined by the terms associated with v„p dLi & Lil 
120061) . This essentially derives from the requirement that the 
proton-alpha velocity difference vector be aligned with the in- 
stantaneous magnetic field. As a consequence, these terms have 
no contribution to the overall angular momentum flux convected 
by the ion flow Lp, which turns out to be small er than in al l 
the models examined in the parameter study bv lLi et aT l (l2007h . 
This, together with the fact that Lp is always positive (i.e., in the 
sense of corotation with the Sun), is at variance with the Helios 
measurements. 

A possible means to reconcile the measurements and the 
model computation is to incorporate the species temperature 
anisotropics. This is because the total pressure tensor P - 2 j Ps 
summed over all species s participates in the problem of angular 
momentum transport via the component - P" - P-*- where 
II and ± are relative to the magnetic field B (see e.g., IWebej 
1 19701 hereafter referred to as W70). While the overall loss rate 
per steradian L may not be significantly altered, the azimuthal 
speed of the solar wind and therefore the particle part of L may 
be when compared with the isotropic case. Note that in the treat- 
ment of W70 the solar wind was seen as a bulk flow and the ion 
species are not d isting uished. On the other hand, the formula- 
tion by iLi & Lil ( |2006|) did not take into account the pressure 
anisotropy, which is a salient feature of the velocity distribu- 
tion functions for both prot ons and alpha particl es as revealed 
by the Helios measurements dMarsch et al.ll982allbb . It therefore 
remains to be seen how introducing the pressure anisotropy in- 
fluences individual ion azimuthal speeds. Moreover, the simple, 
prescribed functional form for assumed in W70 needs to be 
updated in light of the more recent particle measurements. 

The aim of the present paper is to extend the W70 study in 
three ways. First, we shall follow a multicomponent approach 
and examine the angular momentum transport in a solar wind 
comprising protons, alpha particles and electrons where a sub- 
stantial proton-alpha particle velocity difference exists. Second, 
although following W70 we use a prescribed form of for sim- 
plicity, this prescription is based on the Helios measurements, 
and also takes into account other in situ and remote sensing 
measurements. Third, unlike W70 where the model equations 
are restricted to the equatorial plane, the equation set we shall 
derive is appropriate for a rather general axisymmetrical, time- 
independent, multicomponent, thermally anisotropic flow ema- 
nating from a magnetized rotating star. We note that a similar set 
of equat ions, which was also restricted to radial flows, was de- 
rived bv llsenberg| ( ll984 who worked in the corotating frame of 
reference and neglected the azimuthal dynamics altogether. The 
functional dependence on the radial distance and flow speed of 
the magnetic spiral angle was prescribed instead. His approach 
is certainly justifiable for the present Sun, but a self-consistent 
treatment of the azimuthal dynamics is required when flows from 
other stars are examined. This is because many stars either have 
a stronger magnetic field or rotate substantially faster than the 
Sun. 

The paper is organized as follows. We start with section |2] 
where a description is given for the general multifluid, gy- 
rotropic transport equations, based on which the azimuthal dy- 
namics of the multicomponent solar wind is examined. Then sec- 
tion [3] describes the adopted meridional magnetic field and flow 
profiles. The numerical solutions to the angular momentum con- 
servation law are given in section|4] In section|5] we shall discuss 



how examining the angular momentum transport in a multicom- 
ponent solar wind can also shed some light on the spectra of ion 
velocity fluctuations induced by Alfvenic activities. Finally, sec- 
tion|6] summarizes the results. The equations of and a discussion 
on the poloidal dynamics are presented in the appendix. 

2. Mathematical formulation 

Presented in this section is the mathematical development of the 
equations that govern the angular momentum transport in a time- 
independent solar wind which consists of electrons (e), protons 
(p) and alpha particles (a). Each species s {s - e, p, a) is char- 
acterized by its mass m,, electric charge gj, number density «j, 
mass density ps - iisnis, velocity v.,, and partial pressure ten- 
sor p.5. If measured in units of the electron charge e, may be 
expressed by e , - Z^e with s - 1 by definition. 

To simplify the mathematical treatment, a number of as- 
sumptions have been made and are collected as follows: 

1. Symmetry about the magnetic axis is assumed, i.e., d/dcp = 
in a heliocentric spherical coordinate system (r, 9, (p). 

2. The velocity distribution function (VDF) of each species 
is close to a bi-Maxwellian, and the pressure tensor is gy- 
rotropic, i.e., = pf I + (pi - pf)bB, where I is the unit dyad 
and b is the unit vector along the magnetic field B. The tem- 
peratures pertaining to the degrees of freedom parallel and 
perpendicular to B follow from the relation pi'''' - nskgTl:^, 
where kg is the Boltzmann constant. 

3. Quasi-neutrality is assumed, i.e., iie - J^k^k^k ■ 

4. Quasi-zero current is assumed, i.e., Ve = Tjk^knkVk/iie (k = 
p, a), except when the reduced meridional momentum equa- 
tion is derived. 

2.1. Multi-fluid equations 

The equations appropriate for a multi-component solar wind 
plasma with gyrotropic species pressure tensors may be found 
by neglecting the electron ine rtia (nig = 0) in the equations given 
by iBarakat & Schund (Il982l). Following th e same procedure as 
given in the appendix A.I in lLi & Lil(l2006l) . one may find 
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where the subscript s refers to all species (s - e, p, a), while 
k stands for ion species only (k = p, a). The gravitational con- 
stant is denoted by G, Mq is the mass of the Sun, and c is the 
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speed of Ught. The momentum and energy exchange rates due 
to the Coulomb collisions of species s with the remaining ones 
are denoted by 6MJ6t and 6E\'^ jdt, respectively. The third- 
rank tensor Q,, together with the heat flux vectors ^j""" associ- 
ated with parallel and perpendicular degrees of freedom, arises 
from the deviation of species VDFs from an exact bi-Maxwellian 
jBarakat & Schunklll982l) . Moreover, //f'""" stands for the heat- 
ing rates applied to species s in the parallel and perpendicular 
directions from some non-thermal processes. They may be de- 
termined by assuming that the heating de rives from the dissi- 
pation of Alfven-ion cyclotron waves (e.g.. lHollweg & Isenberd 
l2002h. or more simply in some ad hoc fashion such as employed 
in Leer & AxfordI (Il972l) . The operators Vy and are defined 
by V|| -bb-V and Vj^ = V - Vy, respectively. 

In Eq.(|2]i, the subscript j stands for the ion species other than 
k, namely, j - p for k - a and vice versa. As can be seen, in 
addition to the term (V x B) x B, the Lorentz force possesses a 
new term in the form of the cross product of the ion velocity dif- 
ference and magnetic field. Physically, this new term represents 
the mutual gyration of one ion species about the other, the axis 
of gyration being in the direction of the instantaneous magnetic 
field. Furthermore, Equation (|5]l is the time-independent version 
of the magnetic induction law, which states that the magnetic 
field is frozen in the electron fluid. It may be readily shown that 
the effects of the electron pressure gradient, the Hall term, and 
the momentum exchange rates as contained in the generalized 
Ohm's law can be safely neglected given the large spatial scale 
in question (A formal evaluation of t he different terms can be 
found in section 2.1 of lLi et alj ( l2006h ). 

To proceed, we choose a flux tube coordinate system, in 
which the base vectors are {e;, e^j, e^}, where 

ei = Bp/\Bp\, = 60 X e/, 

with the subscript P denoting the poloidal component. 
Moreover, the independent variable I is the arclength along the 
poloidal magnetic field line measured from its footpoint at the 
Sun. This choice permits the decomposition of the magnetic field 
and species velocities as follows. 



B - Biei + B^e^, 



(6) 



where s - e, p, a. From the assumption of azimuthal symme- 
try, and the assumption that the solar wind is time-independent, 
one can see from the poloidal component of equation (|5]) that 
Vep should be strictly in the direction of Bp. In other words, 
VeN = to a good approximation. Now let us consider the 
component of the momentum equation (|2]i. Since the frequen- 
cies associated with the spatial dependence are well below the 
ion gyro-frequency Q.^ = {ZkeBi)/ (111^0) [k - p,a), from an 
order-of-magnitude estimate one can see that Ivj^ - vtml <s: \vii^\. 
Combined with the fact that v^at - 0, this leads to that both Vp^ 
and VaN should be very small and can be safely neglected unless 
they appear alongside the ion gyro-frequency. With this in mind, 
one can find from the component of equation (|2]i that 



B 



— {Val-Vp,). 



(7) 



That is, the ion velocity difference is strictly aligned with the 
magnetic field. This alignment condition further couples one ion 
species to the other 

The fact that v,^ {s - e,p,a) is negligible means that the 
system of vector equations may be decomposed into a force bal- 
ance condition across the poloidal magnetic field and a set of 



transport equations along it. In the present paper, however, we 
simply replace the force balance condition by prescribing an an- 
alytical meridional magnetic field configuration. Moreover, we 
examine in detail only the azimuthal dynamics, leaving a brief 
discussion on the poloidal one in the appendix. 

2.2. Azimuthal dynamics 

The (f) component of the magnetic induction law Q gives 



— {BijiVeP - VetpBp^ 



= 0. 



(8) 



Now that Vep - Veiei, one may readily integrate Eq.® along a 
magnetic line of force to yield 



= AaR + —Vei- 



(9) 



Here R - r sin 6* is a geometrical factor to be evaluated along 
a given line of force (see FiglT]i, and An is a constant of inte- 
gration and should be identified as the angular rotation rate of 
the footpoint of the magnetic flux tube. Taking into account the 
alignment condition (|7]i, one may find that 



= AqR + —Vsi 
Bi 



(10) 



where s - e,p, a. Therefore in a frame of reference that coro- 
tates with the Sun, the velocities of all species are aligned with 
the magnetic field. 

Another equation that enters into the azimuthal dynamics is 
the (p component of the total momentum. In the present case, it 
reads 




PkVkl 



[Rvk^) 



l-—\RB, 



= 0, 



where 

pA ^ pii _ p\\,± ^ p\u^ 



(11) 



(12) 



and the prime ' = e/ ■ V is the directional derivative along the 
poloidal magnetic field. 

For a time-independent flow pkvu /Bi - const. It then follows 

that 



R 



Vp0 + riVa4, - 



BiB^ 

47TppVpl 



1 - 



B'- 



^A, 



(13) 



where the constant 77 - {paVai)liPpVpi) is the ion mass flux ratio, 
and Ai is a constant of integration. Physically, Ai is related to 
the angular momentum loss rate per steradian L by 



L = MpAL, 
where 

B[]7 2 

Mk = PkVki—r^, 
Bi 



(k - p, a) 



(14) 



(15) 



is the ion mass loss rate per steradian scaled to the Earth orbit 
r£ = 1 AU, with Bie denoting the strength of the poloidal mag- 
netic field at r£. It follows that the angular momentum loss rate 
of the Sun due to the solar wind L = AttL if L is independent 
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of colatitude. Equation ( fT3] l shows that L consists of the contri- 
butions due to individual ion angular momenta L^, the magnetic 
stresses Lm and the total pressure anisotropy Lam, where 



[Z^^, Zv^, Lcmi\ 

Bn 



4jt 



B2 



(16) 



with k - p,a. 

Substituting Eq.dTOli into ( fT3b , one may find 

tanOfM^ -(l -/3*cos^ (D)] 
Al 



- e 



(17) 



where tan cD = B^/Bi defines the magnetic azimuthal angle CD, 
and 



Ml ^ Ml + Ml, Ml = 



(k = p, a). 



(1 



Vpl 



(18) 



By definition, Mj is the combined poloidal Alfvenic Mach num- 
ber, which involves both ion species. For a typical solar wind, 
between 1 Rq and 1 AU there exists a point where Mj = 1, 
which is to be called the Alf ven point and d enoted by r„. 

As discussed in detail bv lLi&Lil(l2006h . when species tem- 
perature anisotropy is absent (P^ = and therefore (^^ - 0), for 
Eq.lfTTTi to possess a solution that passes smoothly through ra the 
two constants Ai and An have to be related by 



AL^{l+7j)AnRi, 



(19) 



where the subscript a denotes quantities evaluated at the Alfven 
point. When /S'^ is not zero, a direct relation between Al and 
An is not as obvious since now Eq.dTTli becomes cubic in tan <1>. 
Nevertheless, one may write A/, as A^ - AA^q, where A^o is 
determined through Eq.(fT9]l and therefore A stands for the cor- 
rection due to a finite j6'^. It then follows that 



C3 tan^ C) + C2 tan^ O -i- ci tan O + C2 = 0, 



(20) 



where 



C3 = Mj. - 1,C2 = e 



Given the meridional flow profiles along a prescribed magnetic 
field line, Eq.(|20t possesses only one real root at locations far 
away from r^. However, in the vicinity of r^,, there exists in gen- 
eral three real roots and they diverge near r^. The requirement 
that there exists a unique solution that is smooth from 1 Rq out 
to 1 AU determines A tWeber & Davis..l970;,Weber..l97a) . 

3. Meridional magnetic field and flow profiles 



In principle, one needs to solve Eqs. dA.ll ) to jAAj together with 
Eq.(fT7]i simultaneously to gain a a quantitative insight. In the 
present paper, we refrain from doing so because from previ- 
ous experience it proves diflicult to yield the flow profiles that 
satisfactorily reproduce in situ measurements such as made by 
Helios. Take the proton-alpha speed difference Va,, in the fast so- 
lar wind for example. It is observationally established that Vap 



5F -rrrp-^ 



o 




Fig. 1. Adopted meridional magnetic field configuration in the 
inner corona. Here only a quadrant is shown in which the mag- 
netic axis points upward, and the thick contours labeled F and S 
delineate the lines of force along which the fast and slow solar 
wind solutions are examined, respectively. Also shown is how to 
define the geometrical factor R, and the base vectors e/, ej^ and 
of the flux tube coordinate system (see section|2ll. 



closely tracks the local Alfven spe ed in the heliocentric range 
r > 0.3 AU dMarsch et alj|l982a|) . So fai- this fact still poses 
a theoretical challenge: adjusting the ad hoc heating parame- 
ters, or fine-tuning the cyclotron reson ance mechanism is u nable 
to produce such a behavior (see, e.g. iHu & Habballll999[) . We 
therefore adopt an alternative approach by prescribing the back- 
ground meridional flow profiles that mimic the observations and 
then examining what consequences the species anisotropies have 
on the azimuthal dynamics. 



(21) 3.1. Background meridional magnetic field 



For the meridional magnetic field, we adopt an analytical model 
given bv lBanaszkiewicz et al.l ( Il998h . In the present implemen- 
tation, the model magnetic field consists of the dipole and 
current-sheet components only. A set of parameters M - 2.2265, 
Q - 0, K - 0.9343 and oi - 1.5 are chosen such that the last 
open magnetic field line is anchored at heliocentric colatitude 
6 = 50° on the Sun, while at the Earth orbit, the meridional 
magnetic field strength Bi is 3y and ind ependent of colatitude 6, 
consistent with Ulysses measurements dSmith & Baloghiri995h . 

The background magnetic field configuration is depicted in 
Figdl where the thick contours labeled F and S represent the 
lines of force along which we examine the fast and slow solar 
wind solutions, respectively. Tube F (S), which intersects the 
Earth orbit at 70° (89°) colatitude, originates from 6 - 38.5° 
(49.4°) at the Sun where the meridional magnetic field strength 
B, is 3.93 (3.49) G. 



Bo Li & Xing Li: Angular momentum transport in multicomponent winds 



5 



1000 
^ 800 






(a) : 


600 

g 

^ 400 
^ 200 


- 
- 

// 




- 
■ 








- 








10^ 




r ' 

- - a 


(b) ; 


10^ 














^ 10' 




T V 




10^ 
10' 




.... 1 






_ _ 




1 


10 


100 



3 400 




10^ 



10 

r(Ro) 



100 



Fig. 2. Radial distribution between 1 /?© and 1 AU of the adopted meridional flow parameters for the fast (left column) and 
slow (right) solar wind, (a) and (c), the meridional flow speeds of protons {Vpi) and alpha particles (vq,;). (b) and (d), the ion 

temperatures t\ (dotted lines), (dashed lines), and - {T\ + 2Tt)l3 (solid lines) where k - p, a. The construction of r!''^ is 



described in section [X2l T he error bars in (b) and (d) represent the uncer tainties of the UVC S measurements of the eff'ective proton 
temperature as reported by iKohl et alj d 19981) for a coronal hole, and bv lFrazin et alj ( l2003l) for a streamer, respectively. Note that 
both measurements are typical of solar minimum conditions. Moreover, the asterisks in (a) and (c) denote the Alfven point, where 
the meridional Alfvenic Mach number (defined by Eq.lfTSIl) equals unity. 



3.2. Prescribed meridional flow profiles 

The background meridional flow parameters are found by adopt- 
ing a three-step approach described as follows: 

1 . Using some ad hoc heating parameters, we solve along flux 
tube F (S) the i sotropic version of Eqs. dA.ll i to ( |A.4t (see 
Eqs.(8) to (10) in lLi & Lil ( |2007|) for details) to yield the dis- 
tribution between 1 Rq and 1 AU of the ion number densi- 
ties nit and meridional speeds Vki {k = p, a), as well as the 
isotropic species temperatures Tg (s = e, p, a) for the fast 
(slow) solar wind. Specifically, th e heating rates are of the 
same format as in section 3.2 in iLi & 111 ( 120081) . To gen- 
erate the fast and slow solar wind solutions, the parame- 
ters [Fe (in erg cm"^ s"'), 1^ (in Rq), x] are chosen to be 
[1.9,2.2,2.2] and [1,1.8,3.7], respectively. By simply ad- 
justing the heating parameters it proves difficult to produce a 
reasonable Tp profile in that if Tp in the inner corona is close 
to observations then Tp at 1 AU is usually only a fraction of 
the typically measured values. Moreover, the derived speed 
difference Vapj - Vai - Vpi varies little between 0.3 and 1 AU, 
in contrast to the Helios measurements. Therefore some ad- 
ditional steps are employed to make the flow profiles more 
realistic. Specifically, all the parameters except itp and Vpi are 
required to undergo a smooth transition from the profiles for 
the region r < 0.3 AU derived so far to those specified in 
next step for the outer region. 



2. The desired profiles, given in Table [T] for Vapj, Tp and Ta 
for the region r > 0.3 AU (denoted by subscript o) are based 
on the in situ measurements to be detailed shortly. Once Vapj 
is known, the meridional alpha speed v„/ is given by Vpi + 
Vapj, and the alpha density /Iq, by {na)i{vai)i jvai, where the 
subscript / denotes the values obtained in the first step. The 
distributions of rik, Vki {k - p,a) and Ts {s - e,p,a) thus 
constructed are for the isotropic model. 

3. Now the ion temperatures rl'""" can be constructed by pre- 
scribing the temperature anisotropy = T^l/T^ (k - p,a). 
Note that the electron temperature is assumed to be isotropic. 
For the region within several solar radii, F^ is required to 
decrease with r from 1 at 1 Rq, where the Coulomb self- 
collisions are still frequent enough to suppress a temperature 
anisotropy, to some value less than unity. This inner profile 
is not directly constrained by observations but constructed 
by noting that the processes operational in the inner corona 
tend to heat the ions preferentially in t he perpendicular di- 
rection (e.g.. lHollweg & Isenberg|l2002l) . On the other hand, 
for r > 0.3 AU, F<. follows a power law dependence on r 
with the exponent determined by the Helios measurements 
(see Table[TJ. Specifying the temperature anisotropics of pro- 
tons and alpha particles at 1 AU, YpE and TaE, determines 
Tk in the outer region. The inner and outer profiles are then 
connected smoothly to yield the desired F^. The tempera- 
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tures T^'""" follow from the relations — 3Tk/(2 + Tt) and 



A detailed description of Table [T| is necessary. Note that 
throughout this table x = r/r^ where ve - I AU. Let us first 
focus on the adopted values for the fast solar wind. For the 
isotropic proton and alpha temperatures at 1 AU, we adopted 
the typical values of - 2.8 x 10^ K and Ta - 5Tp (see 
e.g.. lSchwennll990l:lMcComas et al.l2000|) (hereafter Sc h90 and 
McOO). Furthermore, Figures 18 and 19 in lMarsch et ai](ll982bh 
(hereafter M82b) indicate that rj| oc x"°■^^ and oc x"' °^ A 
power law dependence for Tp of Tp oc x ' is therefore consis- 
tent with such a behavior, and also consistent with the Ulysses 
measurements (see Table 2 in McOO). Furthermore, Figure 5 
in iMarsch et all (Il982al) (hereafter M82a) indicates that oc 
x"'-'^, and oc jc"' -'^. A profile of r„ oc x"'-^ is consis- 
tent with this behavior, but differs substantially from that mea- 
sured by Ulysses, which yields that Ta oc (see Table 2 in 
McOO). Moving on to the slow solar wind, we note that values of 

= 5.5 X 10^ K and r„ = 1.7 x 10^ K are typically found at 
1 AU (see e.g., Sch90). In addition. Figures 18 and 19 in M82b 
indicate that Tp oc x"' "^, and Tp oc x""'^. Therefore we adopted 
a Tp profile of Tp oc x"°^'*. On the other hand, we adopted a 
profile for Ta in the form Ta oc x"" which is consistent with 
the measured alpha temperature anisotropics which indicate that 
tI oc x-°•'*^ and T^ oc x-'^^ (see Fig.5 in M82a). 

In this study VpE and FaE will serve as free parameters. 
The Helios measurements indicate that TpE ~ 1.2 + 0.3 and 
FaE ~ 1.3 ± 0.6 for the fast solar wind with Vpi > 600 km s"', 
while X 1.7 + 0.7 andFaf 1.4+0.6 forthe slow solar wind 
with Vpi < 400 km s"' ( iMarsch et al. 1982a,b,). Theoretically, 
one may expect that the [TpE, FaE] pair may not occupy the 
whole rectangle bounded by the given values in the FpE-FaE 
space, since too strong an anisotropy can drive the system un- 
stable with respect to a number of instabilities when the plasma 
yS is comparable to unity. Given that the lower limit of FpE or r^.^ 
is only slightly lower than 1 , the ion-cyclotron instab iUty can be 
show n to be unlikely to occur (see, e.g., Eq.(3) in iGarv et aTl 
Il994l) . However, the firehose instability may be relevant since 
it happens when P" is sufficiently larger than and - 
^nPyB^ > 1. Note that the alpha particles with a non-negligible 
abundance drifting relative to protons may complicate the situa- 
tion considerably given that in addition to the firehose, electro- 
magnetic ion/ion instabilities may also be relevant and the occur- 
rence of such insta bilities is not restricted to the c ases where the 
parallel /3 is large (iHellinger & Travnfcekll20()6l) . Nevertheless, 
we only compare the modeled [rp,rQ,] with the non-resonant 
firehose criterion such as fou nd via the dispersion relation of 
Alfven waves (see Eq.(23) in IIsenberg| |T984). Specializing to 
an electron-proton-alpha plasma, the dispersion relation dictates 
that instability occurs when 1 - > 2(1 - XpX„)/yS" where 

Xk = (pk/p)(vap/vA) (k = p,a) with va = Bj ^[Anp being the 
Alfven speed determined by the bulk mass density p - pp + pa- 
Using this criterion it is found that the modeled flow profiles 
are all stable with the only exception being for the segment 
r > 195 /?0 in the fast wind with the largest values of FpE and 

Figure |2] gives the radial distributions between 1 Rq and 
1 AU of the flow parameters for the fast and slow solar wind 
in the left and right panels, respectively. Figures |2^ and |2}; de- 
pict the meridional ion speeds Vpi and Vai, while the ion tem- 
peratures rj) (the dotted curves), T^ (dashed) and T^ (solid) are 
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given in FigslSJ) and|2]i {k - p, a). The values for the tempera- 
ture anisotropy adopted for the construction are F^ ^ = 1 .5 and 
Fa,E = 1.9 for the fast wind, and Fp_E - 2.42 and Fa.E = 2 for 
the slow wind. In Fig|2j3, the error bars represent the uncertain- 
ties of the UVCS measurements for the proton eff'ecti ve temper- 
ature, made for a polar coronal h ole as reported bv iKohl et al.l 
(Il998l) . Similar measurements bv iFrazin et al.l ( l2003h along the 
edges of an equatorial streamer are given in Figl2}l. Moreover, 
the asterisks in Figs|2h and|2}; mark the location of the Alfven 
point as defined by Eq.lfTSll. 

For the fast (slow) solar wind it is found that at 1 AU the 
meridional proton speed Vpi is 607 (304) km s the proton flux 
iipVpi is 2.8 (3.84) in units of 10*^ cm s"^, the alpha abundance 
Ha/rip is 4.56% (3.6%), and the meridional component of the 
proton-alpha velocity difference Vapj is 23 (5) km s"'. These 
values are c onsistent with in situ m easurements such as made 
by Ulysses jMcComas et al.l |2000|) . Moreover, the fast (slow) 
solar wind reaches the Alfven point at 10.7 (13.3) Rq, beyond 
which Vpi increases only slightly with increasing r. On the other 
hand, for r > 0.3 AU the meridional alpha speed v„/ decreases 
rather than increases with r as a consequence of the prescribed 
Vapj profile. If examining the ratio of Vapj to the meridional 

Alfven speed vai - Bi/ ^jAnp, one may find that for the fast 
solar wind this ratio decreases only slightly from 0.98 at 0.3 AU 
to 0.82 at 1 AU, while for the slow wind it shows a substantial 
variation from 0.88 at 0.3 AU to 0.29 at 1 AU. The modeled 
Vapj I Vai can be see n to agree with the H elios measurements as 
given by Fig. 11 of IMarsch et"aLl ( Il982al) . Note that a value of 
v<yp.i - 49 km s"' at 0.3 AU is not unrealistic for slow solar 
winds, even larger valu es have been found by Helios 2 when 
approaching perihelion jMarsch et alJll98Tl) . Moving on to the 
temperature profiles, one may see that the Tp profiles inside 5 Rq 
are in reasonable agreement with the UVCS line-width measure- 
ments for both the fast and slow solar wind. 



4. Numerical results 

Having described the meridional magnetic field and flow pro- 
files, we may now address the following questions: to what ex- 
tent is the total angular momentum loss of the Sun aff'ected by 
the ion temperature anisotropics? and how is the angular mo- 
mentum budget distributed among particle momenta, the mag- 
netic torque, and the torque due to ion temperature anisotropics? 
To this end, let us first examine the fast and then the slow so- 
lar wind solutions. In the computations, we take An = 2.865 x 
10"^ rad s which corresponds to a sidereal rotation period of 
25.38 days. 
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Fig. 3. Radial distributions of (a) the proton azimuthal speed 
(b) the alpha one \'a^, and (c) various contributions to the 
angular momentum budget in an e - p - a solar wind with ion 
temperature anisotropies. In (a) and (b), the profiles derived for a 
solar wind with identical flow parameters where ion temperature 
anisotropies are neglected are given by dashed lines for com- 
parison. Panel (c) depicts the individual ion angular momentum 
fluxes Lp and L„, their sum Lp, and the fluxes associated with the 
magnetic stresses Lu, and with the temperature anisotropies Lo,„ 
(see Eq.(fT6]l). The dash-dotted lines represent negative values. 



4. 7 . Fast solar wind 

Figure |3] presents the radial profiles of (a) the proton azimuthal 
speed Vp0, (b) the alpha one v„0, and (c) the ion angular mo- 
mentum fluxes Lf: (k = p, a), their sum Lp, the flux due to the 
magnetic torque Lm, and that due to temperature anisotropies 
Lani (see Eq. (fT6] l). Note that the dash-dotted curves in Fig|3}; 
plot negative values. In Figs[3^ andjSj', the ion azimuthal speeds 
for the isotropic model with identical meridional flow parame- 
ters are given by the dashed lines for comparison. The fast wind 
profile corresponds to FpE -1.5 and FaE = 1-9. 



For the chosen TpE and TaE, it is found that A - 1.058. 
Consequently, the total angular momentum loss rate per stera- 
dian L is 1.8 (here and hereafter in units of 10^^ dyne cm sr"') 
in the anisotropic case, and is only modestly enhanced compared 
with the isotropic case, for which L = 1 .7 1 . Furthermore, Figsl3^ 
and[3j5 indicate that the radial dependence of the ion azimuthal 
speed Vp^ or Vq,0 in the anisotropic model is similar to that in the 
isotropic one. For instance, both models yield that with increas- 
ing distance the alpha particles develop an azimuthal speed in 
the direction of counterrotation with the Sun: Vq,^ becomes neg- 
ative beyond 7.95 (8.35) Rq in the anisotropic (isotropic) model. 
The difference between the isotropic and anisotropic cases be- 
comes more prominent at large distances where becomes in- 
creasingly significant, as would be expected from Eq.lfTTli. Take 
the values of Vp^ and v„0 at 1 AU. The isotropic (anisotropic) 
model yields that Vp^ = 2.54 (3.46) km s"' and that v„0 = -12.6 
(-1 1.7) km s"' at 1 AU. Note that the changes introduced to the 
ion azimuthal speeds by pressure anisotropies (~ 0.9 km s"' for 
both protons and alpha particles) play an important role in the 
distribution of the angular momentum budget L among differ- 
ent contributions, as shown by Fig|3};. The proton contribution 
Lp exceeds Lm for r > 57 Rq and Lp attains 5.07 at 1 AU, sig- 
nificantly larger than the magnetic part Lm = 1-48 at the same 
location. In fact, the overall particle contribution Lp, which in- 
creases with distance, overtakes the magnetic contribution Lm 
from 170.5 Rq onwards, despite the fact that the alpha contribu- 
tion tends to offset the proton one. The dominance of Lp over Lm 
happens in conjunction with the increasing importance of Lani, 
the flux due to total pressure anisotropy which is in the direction 
of counterrotation with the Sun. In contrast, without pressure 
anisotropies, at 1 AU it turns out that even though a value of 3.73 
is found for Lp, it is almost cancelled by an L„ of -3.49. The re- 
sulting Lp is thus 0.23, substantially smaller than Lm, which is 
nearly identical to the value found in the anisotropic model. This 
contrast between anisotropic and isotropic cases is understand- 
able since it follows from Eq.(fT3ll that, given that the constant 
Al does not vary much from the isotropic to anisotropic model, 
the change of Lp should be largely offset by that of Lq,,, . 

Figure |4] expands the obtained results by displaying the de- 
pendence on FpE and FaE of (a) the factor A, (b) the proton 
azimuthal speed Vp^ and (c) the alpha one Vq,^ at two different 
distances plotted by the different linestyles indicated in (b), as 
well as (d) the constituents comprising the angular momentum 
flux at 1 AU. In addition to the individual ion contributions Lp 
and La, and the magnetic one Lm, the overall particle contribu- 
tion Lp - Lp + La is also given. Note that -La instead of L„ 
is plotted in Fig|4}J. Moreover, the horizontal bars on the left 
of Figsl^^ and|4}; represent the azimuthal ion speeds derived in 
the isotropic case at the corresponding locations for comparison. 
The open circles correspond to the cases where r„£ = 1.3. It 
turns out that at any given FpE each parameter varies monoton- 
ically from the value with FaE - 0.7, represented by the end of 
the arrow, to the value with FaE = 1.9 given by the arrow head. 
In Figl^^ the arrows have been slightly shifted from one another 
to avoid overlapping. 

From Figll^ one can see that A decreases with increasing 
FpE or FaE, ranging from 1.101 at the upper left to 1.058 at the 
lower right corner. The deviation of A from unity, albeit modest, 
indicates that the changes introduced in the total angular mo- 
mentum loss due to the ion pressure anisotropies are not neg- 
ligible. From Figs|4j3 and |4j; one can see that between 0.3 to 
1 AU, the magnitude of the azimuthal speeds of both species 
decreases with increasing distance. Furthermore, at either 0.3 or 
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Fig. 4. Values of several parameters as a function of TpE, the 
proton temperature anisotropy at 1 AU. (a) The factor A, the de- 
viation of which from unity represents the correction to the total 
angular momentum loss due to the introduction of ion pressure 
anisotropics; (b) and (c) the proton and alpha azimuthal speeds 
Vp0 and at two different heliocentric distances given by differ- 
ent line styles as indicated in (b); (d) various components in the 
angular momentum flux at 1 AU including individual ion contri- 
bution Lp and La, the overall particle contribution Lp - Lp + La, 
as well as the contribution from magnetic stresses Lm- Note that 
—La instead of Lq, is given in (d). The short horizontal bars in 
panels (b) and (c) represent the azimuthal ion speeds derived in 
the isotropic model for comparison. Furthermore, in panel (b) 
the curve corresponding to 1 AU is slightly shifted relative to 
that for 0.3 AU to avoid the two overlapping each other. The 
open circles correspond to the cases where TaE is fixed at 1.3, 
where F^^ is the alpha temperature anisotropy at 1 AU. At a 
given VpE each parameter varies monotonically from the value 
with TaE - 0.7, represented by the end of the arrow, to the value 
with TaE = 1-9 given by the arrow head. The ranges in which 
TpE and TaE vary are determined from the Helios measurements 
(see text for details). 



1 AU, both \'p^ and Vq.^ increase when TpE or F„£ increases. Take 
the values at 1 AU for instance. One can see that Vp^ ranges 
from 2.36 to 3.46 km s"', while Vq,^ varies between -12.8 and 
- 1 1 .7 km s"' . For the majority of the solutions both Vp^ and Vq.^ 
tend to be larger in the algebraic sense than the corresponding 
values in the isotropic model, which yield 2.54 and -12.6 km s"' 
for protons and alpha particles, respectively. However, at 0.3 AU 
Vp^ or Vaip tends to be smaller in the anisotropic than in the 
isotropic case. Now Vp^ and Vq.^ vary in the intervals [2.37, 3.77] 
and [-15.8,-14.2] km s"', respectively. For comparison, the 
isotropic model yields a Vp^ (v^^) of 3.33 (-14.7) km s"'. Now 
let us examine the specific angular momentum fluxes Lp, La and 
Lm at 1 AU. Figure|4}i indicates that Lm has the weakest param- 
eter dependence, which is easily understandable given that to a 
good approximation tan a; -AnR/vu where k may be taken to 
he p or a (see Eq.dTOt). Besides, the parameter dependence of La 
is rather modest, varying by < 10% from -3.55 to -3.23 when 
TpE or TaE changes. On the other hand, Lp changes substantially, 
ranging between 3.46 and 5.07. Hence the overall particle con- 
tribution Lp also shows a significant parameter dependence. In 
particular, Lp may exceed Lm when TpE > 1.3. For the solu- 
tions examined, Lp can be found to be positive and attain its 
maximum of Lp - 1.84 when [TpE, TaE] = [1-5, 1.9]. Only for 
the lowest values of TpE and TaE can one find a negative Lp of 
-0.084. Moreover, the protons always show a partial corotation, 
i.e., Lp > 0. From this we conclude that the ion temperature 
anisotropics are unlikely the cause of the tendency for Lp or Lp 
to be negative for the fa st solar wind as indicated by the Helios 
measurements jPizzo et al...l983;.Marsch & Richter..l984i) . 



4.2. Slow solar wind 

Figure |5]presents, in the same fashion as FiglH the dependence 
on TpE and TaE of various quantities derived for the slow so- 
lar wind. A comparison with Fig|4] indicates that nearly all the 
features in Fig|5] are reminiscent of those obtained for fast so- 
lar wind solutions. However, some quantitative differences exist 
nonetheless. For instance, when TpE is held fixed, all the exam- 
ined parameters for the slow wind vary little even though TaE 
changes considerably from 0.8 to 2. In contrast, the parameters 
for the fast wind show an obvious F^.^ dependence. This differ- 
ence can be largely attributed to the fact that in the slow wind 
the ions are substantially cooler than in the fast wind. Figure|5t 
shows that A ranges from 0.94 to 1.016. In other words, rela- 
tive to the isotropic case, the solar angular momentum loss rate 
per steradian in the anisotropic models may be enhanced or re- 
duced by up to 6%. If examining Figs|5}5 and|5};, one may find 
that at both 0.3 and 1 AU, the azimuthal speeds of both ion 
species, Vp^ and Vq,0, are larger algebraically in the anisotropic 
models than in the isotropic one. The difference between the two 
is more prominent at 0.3 AU, where the isotropic model yields 
that [Vp^,Va^] - [3.49,-18.1] km s"', whereas the anisotropic 
models yield that with increasing TpE, Vp^ increases from 3.76 
to 5.04 km s"', and Vq.^ varies between -17.8 to -16.3 km s"'. 
As for the ion azimuthal speeds at 1 AU, one can see that vary- 
ing TpE leads to a Vp^ varying between 1.18 and 1.72 km s"', and 
a V|j0 ranging from -5.85 to -5.31 km s"'. The coiTcsponding 
changes in the specific ion angular momentum fluxes are shown 
by Fig|5}J, which indicates that the proton one Lp increases with 
increasing TpE from 2.52 to 3.67, and likewise, the alpha one La 
increases from -1.83 to -1.66. On the other hand, the flux as- 
sociated with magnetic stresses Lm hardly varies, and a value of 
3.36 can be quoted for all the models examined. Therefore in the 
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Fig. 5. Similar to Figure |4]but for the slow solar wind. Here the 
open circles correspond to the cases where VaE is fixed at 1 .4, 
and the arrow represents how the specific parameter varies at a 
given FpE when FaE increases from 0.8 to 2.0. 



parameter space explored, Lm may be smaller than Lp, which 
is however offset by the alpha contribution that is always in the 
direction of counter-rotation to the Sun. In fact, the alpha contri- 
bution is so significant that the overall particle contribution Lp 
never exceeds Lm- In other words, incorporating ion temperature 
anisotropy cannot resolve the outstanding discrepancy between 
previous models and observations concerning the relative im- 
portance of particle and magnetic contributions in the angular 
momentum budget of the solar wind. 

5. Discussion 



As demonstrated by iLi & Lil (l2008h . the discussion on the an- 
gular momentum transport also allows us to say a few words 
on the frequency spectra S k{f) {k = p,a) of the ion veloc- 
ity fluctuations during Alfvenic activities in the fast solar wind 
in the super- Alfvenic portion where » 1. This is due to 
the well-known change of the properties of Alfvenic fluctua- 
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Fig. 6. Radial dependence of the ratio of the alpha to the pro- 
ton velocity fluctuation amplitudes 5ual6up induced by Alfvenic 
activities in super-Alfvenic portions of the fast solar wind. The 
dashed curves correspond to the isotropic model, while the 
hatched areas give the possible range Sua/Sup may occupy when 
the parameters FpE and FaE vary in the ranges given in text. Both 
the zero-frequency (upper portion) and WKB (lower) estimates 
are given. 



of center of mass evaluated at the Alfv en point r,, (see e.g., 
Heinema nn & Olbertlll980l:lLi & Lil l2008'). For typical fast wind 
parameters, fc ^ 0.5 - 1 x 10"^ s"'. While the fluctuations with 
frequencies f ^ fc are genuinely wave-like and may be de- 
scribed by the WKB limit given the slow spatial variation of 
flow parameters in the region in question, those with / < f- 
behave in a quasi-static manner and may be described by the so- 
lutions to the angular momentum conservation law wh ich also 
governs the zero-frequency fluctuations. As shown by iLi & Lil 
(I2008h who neglected the species temperature anisotropy, in the 
region r > 0.2 AU which will be explored by the Solar Orbiter 
and Solar Probe, the ratio of the alpha to proton velocity fluc- 
tuation amplitude 6ua/Sup can be an order-of-magnitude larger 
for / < fc than for / > fc- Hence one may expect that, if the 
proton velocity fluctuation spectrum 5',,(/) is somehow smooth 
around fc, then the alpha one S aif) will show an apparent spec- 
tral break. Now let us revisit this problem in light of the discus- 
sion presented in this paper and see what changes the pressure 
anisotropics may introduce. 

Restrict ourselves to either the high-latitude region or the 
region inside say 100 Rq such that the magnetic field may be 
seen as radial. Furthermore, suppose that the waves are propa- 
gating parallel to the magnetic field in the empirical fast wind 
profiles detailed in section[3] Figure |6]presents the radial depen- 
dence of dUa/Sup in the region between 40 and 100 Rq for both 
the zero-frequency (upper part) and WKB (lower part) solutions. 
For comparison, the dashed curves represent the corresponding 
results in the isotropic model. To construct Fig|6] all the pos- 
sible values of FpE and FaE have been examined. As a result, 
at any radial location the ratio SudSup varies from model to 
model, and the range in which this ratio may occupy is given by 
the hatched area. The zero-frequency solutions are obtained by 
solving Eq.(l20li. while for hydromagnetic WKB Alfven waves it 
is wefl known that 6ua/6up = |(vp/, - Va) Kvoh - Vp)l, where Vph is 
the wave phase speed and given by (e.g.. lBarnes & Suffblkll97lT: 
llsenberd 19841) 



tions around some fc 
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in which Vcm — PpVp + PaVa is the speed of center of mass, and 
Pk - Pkl P {k - p, a) defines the fractional ion mass density. 

From Fig |6] one can see that the zero-frequency and WKB 
solutions are well separated from each other, in the isotropic and 
anisotropic cases alike. For the isotropic model, 6ua/Sup in the 
zero-frequency case increases monotonically from 3.79 at 40 Rq 
to 4.68 at 100 Rq. On the other hand, in the WKB case it de- 
creases first from 0.22 at 40 Rq and attains its minimum of 0.083 
at 60.3 and then increases to 0.15 at 100 Rq. The difference 
in dUa/Sup between the zero-frequency and WKB solutions may 
be slightly smaller in the anisotropic than in the isotropic case 
for some combinations of [rp£,r„£], but the difference is still 
quite significant. From this we can conclude that, with realistic 
ion temperature anisotropics included, the alpha velocity fluctu- 
ation spectrum S„{f) during Alfvenic activities will also show 
an apparent break near fc, if the proton one S p{f) is smooth 
there. This break is entirely a linear property, and has nothing 
to do with the nonlinearities that may also shape the fluctuation 
spectra. 

6. Summary 

This study has been motivated by the apparent lack of an anal- 
ysis on the angular momentum transport in a multicomponent 
solar or stellar wind with differentially flowing ions and species 
temperature anisotropy. Moreover, there has been an outstand- 
ing discrepancy between available measurements and models 
concerning the relative importance of the particle Lp and mag- 
netic contribution Lm to the solar angular momentum loss rate 
per steradian L. The Helios measurements indicate that for fast 
(slow) solar wind with Vp > 600 (< 400) km s"', tends to be 
negative (positive), with the positive sign denoting the direction 
of corotation with the Sun. Furthermore, Lp tends to be larger 
than Lm in the slow wind. The behavior of Lp derives from that 
of individual ion angular momentum fluxes, Lp and L„, thereby 
calling for a multifluid approach. 

Starting with a general set of multifluid transport equa- 
tions with gyrotropic species pressure tensors, we have de- 
rived the equations for both the angular momentum conserva- 
tion (Eqs.lfTOt and ( l20b in section|2|, and the energy and linear 
momentum balance (Eqs. dA.ft to ( IA.4l i in the appendix). These 
equations are not restricted to radial outflows in the equatorial 
plane, instead they are valid for arbitrary axisymmetrical winds 
that include two major ion species, and therefore are expected to 
find applications in general outflows from late-type stars. To fo- 
cus on the problem of angular momentum transport, we refrained 
from solving the full set of equations governing the meridional 
dynamics. Rather, we constructed, largely based on the available 
in situ measurements, the empirical profiles for the meridional 
magnetic field and flow parameters. Only the ion temperature 
anisotropics are considered, i.e., the electron temperature is seen 
as isotropic. For both the fast and slow solar wind profiles, we 
solved the angular momentum conservation law (Eqs.lfTOli and 
( |20l )) to examine how the azimuthal speeds of protons Vp^ and 
alpha particles v^^, as well as the individual components in the 
solar angular momentum budget are influenced by the ion tem- 
perature anisotropics. To this end, solutions to the isotropic ver- 
sion are obtained for comparison. 

Our main conclusions are: 

1 . From the derived equations governing the energy transport, 
a simple analysis given in the appendix yields that the adi- 
abatic cooling may be considerably influenced with the in- 
troduction of the azimuthal components. Such an influence 



is understandably more prominent in the low-latitude re- 
gions. This means, when modeling the species temperature 
anisotropy, for a quantitative comparison of model compu- 
tations to be made with the near-ecliptic measurements such 
as made by Helios, the spiral magnetic field has to be taken 
into account. 

2. In agr eement with the single-fluid case ( IWeber & DavisI 
ri970t IWebed Il970l) . mcorporatmg species temperature 
anisotropy leads to a situation where the total angular mo- 
mentum loss rate per steradian L is determined by the be- 
havior of the solution to the angular momentum conservation 
law in the vicinity of the Alfven point where the combined 
Alfvenic Mach number Mp = 1 . However, Mp has to take 
into account the contribution from both ion species, as de- 
fined by Eq.®. 

3. Relative to the isotropic case, the introduced species temper- 
ature anisotropy may enhance or decrease L by up to 10%, 
and introduce an absolute change of up to ~ 1 .8 km s ' in in- 
dividual ion azimuthal speeds in the region between 0.3 and 
1 AU. While these changes seem modest, the correspond- 
ing changes in the angular momentum fluxes convected by 
protons Lp or alpha particles L,, may change substantially. 
In contrast, the flux associated with magnetic stresses Lm 
hardly varies. 

4. However, introducing ion temperature anisotropics cannot 
resolve the discrepancy between in situ measurements and 
models. For the fast wind solutions, while in extreme cases 
Lp may become negative Lp always stays positive. On the 
other hand, for the slow solar wind solutions examined, Lp 
never exceeds Lm even though Lm may be smaller than the 
individual ion contribution. This is because, for both the slow 
and fast wind solutions, Lp and L„ always have opposite 
signs. 

5. The discussion on the angular momentum transport has 
some bearing on the ion velocity fluctuation spectra S tif) 
(k = p,a) during Alfvenic activities in the super- Alfvenic 
regions, which are likely to be explored by future mis- 
sions suchas Solar Orbiter and Solar Probe. In agreement 
with|LL&L3(II008) where species temperature anisotropics 
are neglected, an analysis based on the WKB and zero- 
frequency solutions yields that Saif) will show an apparent 
break around some critical frequency fc if S p{f) is smooth 
there. This f- ~ 0.5 - 1 x 10"^ s"' is the well-known fre- 
quency that separates the genuinely wave-like fluctuations 
from quasi-static ones. 
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Appendix A: Derivation of equations governing tlie 
meridional dynamics 

In section |2] we have demonstrated that the vector equations 
governing a time-independent multicomponent solar wind with 
species temperature anisotropy are allowed to be decomposed 
into a force balance condition across the poloidal magnetic field 
and a set of transport equations along it. The azimuthal dynam- 
ics has been discussed in the text, whereas this appendix pro- 
vides some discussion on the poloidal dynamics. In particular, 
we shall derive the equations governing the poloidal motion v*/ 
of ion species {k - p,a), and the species temperatures T^'^ 
{s - e,p,a)\n rather general situations. 

Due to the presence of VkN in the I component of the ion mo- 
mentum equation (O, one may expect that the A^-component of 
Eq.© has to be solved. In fact, there is no need to do so because 
VkN appears only in the difference Vyw - VkN, which may be found 
from the component of Eq.©. Substituting Vj^ — vun into the / 
component of Eq.(|2| will then eliminate the cu mbersome and 
VkN- N ote that this technique, first devised by iMcKenzie et alj 
(II979I) . ensures the con servation of no t only total momentum 
but also total energy (see lLi &Xill2006l) . Specifically, the result- 
ing equations for the poloidal dynamics are 

Vkiivkd' - v^^(ln/?)' + tanO^CTJvi^)' 

- (Q, + tan OQ^) + ^^(In r)' 
r 

+ — {(p!)'-p^[ln(B;sec(l>)]'} 
+ ^{(p[!)'-/5^[ln(B,sec(I))]'U0, 



(A.l) 



(A.2) 



(r|)'+2v,,rf [ln(v,,sec(5)]' 



1 



6E" 
ot 



nsks 

v./(r,^)'-v,,r,^[ln(B,sec<D)]' 



= 0, 



1 



nskB 



, Qs- - SEf 
V.qf + ^:V{bb)-^-Hf 



= 0, 



(A.3) 



(A.4) 



where p^ - pi - pf (s - e,p,a), Ck - 6Mk/6t + 
{Zknklne)5MelSt results from Coulomb frictions, and the prime 
' represents the derivative with respect to the arclength /. When 
deriving Eqs.( |A.3t and ( |A.4t , we have used the fact that the ex- 
pression 

V|| ■ Vj = cos^ ^v\i + sin^ Ovj/(ln/?)' 
+ cos O sin (I)/?(Vj0 //?)', 

may be simplified by expressing Vj^ via the alignment condi- 
tion ( [Tol l, the result being 

V|| ■ Vj = Vj/ [ln(vj/secO)]' . 

Similarly, one may find that 

Vi-v, = -v,;[ln(B/secO)]'. 

No w let us compare these equations with those in llsenberg) 
( 119841) . It is straightforward to show that Eqs. ( |A.3t a nd ( |A.4| i are 
equivalent to Eqs.(A3) and (A4) in Ilsenberddl984l) by special- 
izing to a spherically symmetric solar wind and by noting that 



(InB/)' = -2/r. On the other hand, using the alignment condi- 
tion ([Tol l one may find that 



VH (vh)'-v^^ (In /?)'-!- tan O 



[Rvk,p) 



^sec^O 






U J 




[ 2 J 



(A.5) 



Again specializing to a spherical solar wind, one then finds that 
Eq.( IA.2b is equivalent to (A2) in Jsenberg (1984|). It should be 
stressed that although working in a fra me of reference corotat- 
ing with the Sun, as did llsenberd (|I984|) . substantially simplifies 
the algebra, it does not offer the information on the specific form 
of the spiral angle O, whose functional depend ence on the flow 
speeds has to be assumed a priori. In practice, llsenberd d 1 9841) 
assumed that the velocity of center of mass Van is radial in an in- 
ertial frame beyond 10 Rq, which is certainly a good assumption 
for the present slow-rotating Sun. However, from our discussion 
on the azimuthal dynamics, there is in general no guarantee that 
Vcm is radial, and the deviation may be substantial for winds that 
flow from a faster rotating star 

Introducing azimuthal components may influence the ion 
flow speeds Vki both directly and indirectly. The direct conse- 
quence is that azimuthal components may introduce into the re- 
duced meridional momentum equation (IA.2b an effective force 
(see the first three terms). Note that in a corotating frame the 
magnitude of the ion velocity becomes Vki sec <1), from rela- 
tion ( IA.5b one may see that in such a frame all particles move 
in the same centrifugal potential A^R^/2. Therefore in effect the 
introduced force tends to reduce the magnitude of the ion speed 
difference with increasing distance as sec O tends t o incr ease. 
This eff e ct has been explored in detail in iLi & Lil ( |2006|) and 
iLi et al. 1 (I2007h . where it is shown that the influence may play 
an important part in the force balance for the solar wind. In fact, 
introducing solar rotation alone is able to reproduce the Vap pro- 
file measured by Ulysses beyond 2 AU if a proper value of Vap 
is imposed there. On the other hand, Vki may be altered indi- 
rectly by the modified pressure gradient force due to changes in 
the temperatures, which in turn are caused by the changes in the 
heat fluxes (the third term in Eqs. (IA.3b and (IA.4b ) and through 
the adiabatic cooling (the second term). A detailed discussion on 
the former requires a specific form for the heat flux, which is be- 
yond the scope of the present paper As a consequence, we shall 
focus on the latter instead. 

Neglecting the terms in the second pair of square parenthe- 
ses, Eqs.( |A.3b and ( |A.4b give 



oc cos^ O/v^,, oc Bi sec O. 



(A.6) 



Note that the relation governing simply reflects the conserva- 
tion of magnetic moment. Now that in the region say r > 10 Rq 
sec <S> is significant and increases with r, (T^) may be sub- 
stantially reduced (enhanced) relative to the case where O s 0. 
This effect is particularly significant in the near-ecliptic region 
and for the slow solar wind. For instance, restrict ourselves to the 
equatorial plane and consider the region between say 10 Rq and 
1 AU. Suppose x'si remains constant and v.,/ ~ AqRe - 430 km/s. 
Now that roughly speaking tanO x -Anr/x'si, when the spiral 
field is considered, (Tf) at 1 AU is 1 /2 ( V2) times the value 
for a purely radial magnetic field. This suggests that for making 
any quantitative comparison of the modeled species temperature 
anisotropy with the near-ecliptic measurements such as made by 
Helios, the spiral magnetic field has to be considered. 

For completeness, we note that the force balance condition 
across the direction comes from the component of the total 
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momentum, which reads 



^ \Rc '"''dN ) dN\ Sn) 

1 / A7tP^\(B] , d \ 
^{l-^r^W^-Bl—MR 



k 

where Rc - e^ ■ (ei ■ Ve/) is the signed curvature radius of the 
poloidal magnetic line of force. Obviously, this force balance 
condition determines the poloidal magnetic field configuration 
in response to the electric currents associated with the flow. 
This equation, combined with the transport equations along the 
meridional magnetic line of force, may be solved alternately to 
find a self-cons istent solution to the vect or equations by u sing 
the approach bv lPneuman &KoDd (1197 Ih or lSakurail (Il985h . 



